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I. MOTIVATION 

The use of sum rules in quantum mechanics is well 
established and abundantly applied for many different 
systems [l|, 0] . The prominent examples are the transi- 
tions from a given quantum state induced by an electro- 
magnetic multipole operator. For any multipole operator 
acting on an initial state, the sum of all the related transi- 
tion probabilities multiplied by powers of the excitation 
energy are completely determined by the properties of 
the initial state [3, Q • 

The sum rules exist in general for any many-body 
quantum system. Of specific interest are those systems 
where the constituents clusterize, such that the degrees 
of freedom can be divided into the internal ones corre- 
sponding to each cluster, and those associated with the 
relative motion of the clusters Then the different 
multipole operators can be decomposed into terms de- 
pending on the intrinsic coordinates of each cluster, and 
an additional term depending only on the relative coor- 
dinates of the centers of mass of the clusters. This opera- 
tor structure then leaves two sum rules showing the same 
decomposition, the sum rules associated with each indi- 
vidual cluster (depending only on the properties of the 
initial cluster state), plus the cluster sum rule (depending 
on the properties of the few-body initial wave function). 
Examples are found in [6, 7], where the dipole non-energy 
weighted and dipole energy-weighted sum rules are ob- 
tained for many-body systems clusterizing into a two- 
body system. 

When a clusterized system can be properly described 
as a few-body system where the internal cluster degrees 
of freedom are frozen, only the cluster sum rules re- 
main, corresponding to the much smaller Hilbert space of 
ground and excited states of the relative cluster motion. 
This kind of few-body descriptions have been extensively 
used in nuclear physics during the latest 10-15 years in 
connection with halos and weakly bound states in general 



[5|. The most interesting and frequently investigated of 
these systems are approximated by a three-body struc- 
ture. Extensions to excited three-body continuum states 
are now being pursued and attracting a lot of attention 
d, [l3, [m ■ To get accurate three-body wave func- 
tions the Faddeev decomposition with different Jacobi co- 
ordinates is employed in coordinate space computations 
[Tsj . The unavoidable transformation from one set of Ja- 
cobi coordinates to another complicates the structure of 
the cluster sum rules, especially when more than one of 
the three particles is charged. 

The purpose of this work is to generalize the dipole 
two-body cluster sum rule as in Q to three-body sys- 
tems for any multipolarity. Advanced three-body calcu- 
lations employ partial wave dependent pair interactions 
and state dependent three-body potentials, and these 
complications must therefore also be taken into account 
in derivations of the cluster sum rules. Typically, the two- 
body interactions are adjusted independently for each 
partial wave in order to reproduce the available prop- 
erties of the corresponding two-body system, e.g. bound 
state and resonance energies, and phase shifts [3, EBl- 
These interactions are then essentially non-local through 
their angular momentum dependence. Also, it is common 
to introduce effective angular momentum and parity de- 
pendent three-body forces for fine-tuning the crucial to- 
tal energies of the three-body states. Since the multipole 
transition operators carry both angular momentum and 
parity they do not commute, in general, with the angular 
momentum dependent two- and three-body potentials. 
Thus, the energy weighted cluster sum rule should then 
be rederived including the corresponding contributions in 
addition to the usual kinetic energy term. In all cases we 
must include contributions from the three Faddeev com- 
ponents which are expressed in their respective Jacobi 
coordinates. 

In section|lI]we briefly introduce the coordinates used, 
and summarize some important relations and definitions. 
In section lnll we derive the non-energy weighted sum rule. 
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The energy- weighted sum rule is obtained in section [TV] 
which is divided into three subsections corresponding to 
the contributions from the kinetic energy operator, the 
partial wave dependent two-body potentials, and the (to- 
tal angular momentum) J'^-dependent three-body forces. 
As an illustration, in section |V] we investigate the dipole 
excitations in ^He and compare with the available exper- 
imental data. We close the paper with a short summary 
and the conclusions. A few intermediate expressions ob- 
tained in the derivations have been collected in the ap- 
pendix. 



II. THE TRANSITION PROBABILITY 

We assume three clusters with masses rrii and charges 
Zi (j=l,2,3), described by coordinates r^, and with the 
three-body center-of-mass at R. The three sets of mass- 
scaled Jacobi coordinates are {a;,, y^} = {p, a^, fij;. , fij,. }, 
where p is the hyperradius and {a^, flxi, ^yi} are the five 
hyperan gles corresponding to the Jacobi set i, see for 
instance [13[ . The connection between the Cartesian and 
the mass-scaled Jacobi coordinates is given by 



rk 



R 



mjiTik 



'niN (rrij + nik) 



M 



(1) 



(2) 



where M—mi+jrij+mk and tojv is the normalization 
mass. The transformation between different sets of Ja- 
cobi coordinates are given by 



-Xk sm ip^k ~ Vk cos ipik 



(3) 



tanipik = (-l)P 



niiM 



rriimk 



(if i ^ k), and ipn = tt, (4) 



which formally amounts to a rotation depending on the 
mass ratios and the parity (—1)'' of the permutation p of 
{i,j,k}. Eqs.(I3]) and ^ lead to an important relation 
between harmonic polynomials in different Jacobi coor- 
dinates, i.e. 



y^Y^A^yJ = ^(-l)^a;^^(sin(^,fc)^-^y^(cos<^,fe)^ 

t.=o 



Let us consider the initial three-body state \noJoMo), 
where Jq is the total angular momentum with projection 
Afo- All the other needed quantum numbers are collected 
into riQ. The excited states {\nJM)} can be populated 
from the ground state by the electric multipolc operator 



0^ = Y.z^\n-R\^Yx,f,{n 



(6) 



where i runs over the three clusters, or equivalently, over 
the three sets of Jacobi coordinates. 

The transition probability corresponding to this elec- 
tric multipolc operator is proportional to the B{EX)- 
value, i.e. 



B{EX,noJQ — > nJ) 



\{nJM\Ol\nMIo)? 



(7) 



from which the A-multipolc strengths are defined as 



Sm = Y.^E^-' - EorBiEX, no Jo n J), (8) 

nJ 

where Eq is the energy of the initial state, and Enj is 
the energy of the excited state with angular momentum 
J and additional quantum numbers n. 

The values of these multipolc strengths, depending 
only on the properties of the initial state, are known 
as the sum rules. In this work we are concentrating on 
the sum rules with m=0,l, also denoted as non-energy 
weighted and energy-weighted sum rules, respectively. 

For the dipole case (A=l), after inclusion of Eq.© into 
©, and using Eqs.® and (g]), one can see that for three 
particles having equal value of the ratio Zi /rrii the dipole 
operator is zero. This means that for the particular case 
of three identical particles all the dipole strengths Sm in 
Eq.® are zero. 



III. THE NON-ENERGY WEIGHTED SUM 
RULE 

The sum over all transitions can be rewritten pro- 
vided the intermediate set of quantum numbers gives a 
complete description of the (bound and continuum) final 
states, i.e. J2njAi \nJM){nJM\ — 1. We then get: 



S'o = B{EX, no Jo — ^ n J) 



nJ 



J2 MMo I Of,^ \nJM) {nJM\0^\noJoMo) 

nJ fj.AI 

J2MMo\O^JO^^\noJoMo) , 



(9) 



which for a given multipolc operator is entirely deter- 
mined by the properties of the ground state jnoJo-Mo). 

The definition in Eq.®, together with ([2|) and per- 
mits expressing the operator O^^O^ in terms of a single 
set of Jacobi coordinates, leading to 



3 A A 



i,k=l 1=0 fj=-A 



{2i + iy.{2x - 2e + ly. 



X-£ X 

1712 — M 



noJoMo\xl '''y^^'''Yx^i,{ny^)Yt,„,-^{Qy^)Yx-e.m2{^Xk)\'noJoMo) 
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TABLE L Non-energy weighted sum rule (So) and the contri- 
bution of the kinetic energy operator to the energy weighted 
sum rule (Sj"^') for a system of three particles with equal 
mass (m) for A=l. The first column gives the number A'' of 
charged particles each with charge z. The symbol () denotes 
expectation value in the initial state, and is the coordinate 
for one of the charged particles. 





So 




1 




|rp - J?p)/(47r) 




2 




\v,-Rf)l{\-.) 
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where the constants 



/mjv (raj -|- toj.) 



M 



(11) 



IV. THE ENERGY- WEIGHTED SUM RULE 

The energy weighted sum rule is most easily obtained 
by evaluating the expectation value of the double com- 
mutator in the ground state, i.e. 

Si = i^(noJoMo|[[0,^t^i/],0^]|noJoAfo> = 
Y,{n^J^M^\0]} HOl- E^O';}Ol\n^J^M^) , (14) 



where is the ground state energy. 

This expression is obtained by inserting the identity 
operator 1 = X^njA/ between H and 
and between and , where {|nJM)} are the com- 
plete set of eigenstates of H with the corresponding set of 
eigenvalues \Enj\- In this way we immediately recover 
the standard definition in Eq.®: 



arise when inserting ([2]) into the definition ([6]). 

The summation over the indexes /i and mi can be made 
analytically [IB], leading to the final expression for the 
non-energy weighted sum rule: 



( \ \-l l\ 2A + 1 



(12) 



^(noJoMo|:r^-^2/^VA_,,™,(r!,jy;_,,„,(r!,J|noJoMo) 

mi 

The sum {So)diag of the diagonal terms {i—k) in Eq. (fT2|) 
is obtained taking (pu—TT and £=X (which reduce Eqs.Q 
and ([5]) to identities) which leads to: 

iSo),ua, ^^^J2''(JoMo\\r,-RnJoMo) . (13) 

1—1 



For a three-body system containing only one charged 
particle the non-energy weighted sum rule reduces to one 
of the three diagonal terms in Eq. (fT3|) . When more than 
one charged particle enter in the three-body system, the 
full expression ([12]), that contains interferences between 
charged particles, must be used. 

The relevance of the non-diagonal terms can be easily 
seen for a system containing three identical particles with 
mass m and charge z for A=l. In this case the sum of the 
diagonal contributions 3^{\rp — given by Eq. fT^ 
is fully canceled by the non-diagonal terms, such that 
So=0, as expected for three identical particles. 

When only two of the particles with mass m each have 
the charge z, one of the diagonal terms is canceled out by 
the non-diagonal one, and we get Sq = ^{\rp — ilp), 
which is identical to the result when only one particle 
is charged. These results are summarized in the second 
column in table H] 



^1 



nj 



nJ 



Eo)J2\{nJM\0^ 



noJoMo)p. (15) 



According to Eq. (fH)) the energy- weighted sum rule de- 
pends on the multipolc operator, the initial state prop- 
erties, and the hamiltonian. This hamiltonian can have 
a complicated angular momentum dependence of both 
two- and three-body interactions, whose contributions to 
Si in general do not vanish. In particular we shall assume 
two-body interactions that depend on the relative partial 
wave between the two particles, and three-body poten- 
tials depending on the total angular momentum and par- 
ity of the three-body state. 

In the following subsections we evaluate the expression 
(fn)) separately for the different terms of the hamiltonian, 
i.e. the traditional contribution from the kinetic energy 
operator, and the new terms arising from the partial wave 
dependent two-body potentials, and J'^-dependent three- 
body potentials. 



A. Kinetic energy operator 

The kinetic energy operator can be expressed in terms 
of any of the three sets of Jacobi coordinates as 



T 



(16) 



where the two Laplace operators and Aj,j. are asso- 
ciated to the Jacobi coordinates Xk and y^,. Since the 
multipole operator ([6]) only depends on y-coordinates we 
can quickly find the commutator [O^^, T] (Eq. (|Aip ). from 
which we get: 



(17) 



V,, (|r, - R\^Y*jny,)) ■ V,, (|r, - fi|^n,^(f7,J) . 
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Eq.([5]) permits rewriting of Eq. (jl7[) in terms of a single 
set of Jacobi coordinates, leading to: 



niM ^ — ' ^ — ' — ' 

i,k=l f=0 mn 



:\/2A + 1 



\-£ A 




47r(2A + 1)! 



(2£+ 1)!(2A- 2£+ 1)! 



(18) 



The scalar product can now be performed by use of the 
Gradient Formula (jA2p . and after writing the two spher- 
ical harmonics in terms of a single one, and perform- 
ing analytically the summations over angular momen- 
tum projection quantum numbers (details are given in 
the appendix), one gets the following final expression for 
the contribution of the kinetic energy operator to the 



energy- weighted sum rule 



SI 



(T) 



2m]' 



{ciY {ckY ZiZk^^isinip^kY ' {cos (pik)(l9) 
i,k = l e=i 



A^(2£-1)(2A + 1)!(2A-H)« / A-1 £-1 \-l 
(2£)!(2A-2£+l)(2A-2£-hl)! \^ 

'^{noJoMo\xl~'^yl+^~'^Yx-tr,{^^k)Yx-i,n{^yk)\noJaMo) 



where the sum of the diagonal parts (i—k) become 
ft2 A(2A + 1)2 



\ / aiaq 



2m N Att 



(20) 



Y,{c,)^zf{noJoMo\\r, - R\^^-^\noJoMo) 



^26 = J2'i=i ^2b where the index i runs over all the three 
sets of Jacobi coordinates and V^y is the two-body op- 
erator describing the interaction between particles j and 
k. This two-body operator is formally written as: 



Y 



2b 



ix^)P^ 



\^Xi T '^Xi T j Xi ^ '^X i} {^Xi T '^Xi T j Xi ^ ^Xi \ ; (21) 



where Pi represents any spin operator that could enter 
in the two-body potentials. 

The contribution of the full two-body potential opera- 
tor V26 to the second sum rule has then three contribu- 
tions, each corresponding to one of the three two-body in- 
teractions. According to Eq. p^ . the contribution S'J^^'*'' 
(i=l,2,3) from each of them is given by: 

S?"'"^ = E [{noJoMo\0';V^:^0'^\noJoMo) 

-(noJoAfo|0;i^O;il>2bVoJoMo)l , (22) 



where |no JqMo) represents the initial state with total an- 
gular momentum Jq and projection Mq. The quantum 
number hq refers to all other additional quantum num- 
bers necessary to specify this state. 

For each two-body interaction it is now conve- 
nient to write the corresponding ground state wave func- 
tion 5" in terms of the Jacobi coordinates {xi,yj} = 
{p, ai, flxi, ^yi}, and expand it in terms of a set of func- 
tions y:^i>^'''>{n,) 



(23) 



where the constants Ci are given by Eg. ljlip . 

For A=l the expression in Eq. p^ is independent of 
the properties of the initial state. In particular, for three 
identical particles with mass m and charge z the total 
value of S{ is zero, which confirms the result antici- 
pated at the end of section |TT1 When one of these three 
particles has no charge, S\ takes a constant value which 
is the same as the one obtained when only one of the 
three particles with mass m is charged. The precise ex- 

(T) 

pressions of S{ for these particular cases are given in 
the last column of table HI 



B. Partial-wave dependent two-body potentials 

Typically, the two-body interactions are adjusted sep- 
arately for the individual partial waves in order to repro- 
duce the known experimental data for the two-body sys- 
tems. This procedure leads to two-body interactions de- 
pending on the two-body quantum numbers {(-xiSxtJx}- 
The full two-body potential operator takes the form 



where 



6^ " Y''''' 

(24) 

with7i = {K,(.x,,Sxi,jxi,tyi,jyi}-, and with (/);^"' ' {ui) 
being the usual function of the hyperangle Ui entering 
in the definition of the hyperspherical harmonics p^ . 
The functions ([M)) reduce to the usual hyperspherical 
harmonics for particles without spin. 

With the definition (I2ip . the two-body potential opera- 
tor acting on a term of the basis 3^J^f ^'^"(^1;) (written 
in the Jacobi set i) leads to 

v^^y'f'°m = ^/'^"'^"'^'^(x.)p,3^^^''^''(f^o (25) 

Eq (22]) permits to write Eq.fg!]) as: 



{2b,i) 



E 



/^^EE^- 

7. 7,' 



(26) 
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Inserting in the first and second matrix elements the 
unity operator 

7f J"M" 

between O^^ and Fjfc'*' ^'^'^ '^/i^ respectively, 
and making use of Eq. (|25p , we immediately get the final 
expression: 

7^ J"Af" 

When the 1^-functions are independent of the partial 
wave, the equation above can be written in a more com- 
pact way as: 

^r"^ = ^E [dpY^Y^K-'^ipKi'^ip) (29) 

which is trivially zero for the central part of the two-body 
potential (Pi = 1), and for the spin-spin term (O^ does 
not depend on the spin and therefore commutes with the 
spin-spin operator). The same happens for the tensor 
operator, which depends only on coordinates and spin 
operators. For the spin-orbit term {Ix ■ ~ ^+s_ + 
£-S+ + izSz) one has the same result, since or 
iz apphed on is proportional to either O^+i, O^-i, 
or O^, and therefore each of the three terms in ■ 
double commutes with Of^. 

Thus, for two-body interactions independent of the 
partial waves and containing the usual spin operators one 
has ) = 0. 

If the the two-body potentials are partial wave depen- 
dent, but contain only central, spin-spin {Pi=Sj ■ Sk) and 
spin-orbit {Pi—lxi ■ SajJ terms, since these operators are 
diagonal in the basis {{ixi, Sxi, jxn'm'Xi)}, one then has: 

where /^'°' ^ =1 for the central part of the potential, 

= {sxiisx, + 1) - Sj{sj + 1) - Sk{sk + l))/2 for 

the spin-spin part, and fj^^s^. = UxiU^i + 1) + 
1) — Sxiisxi + l))/2 for the spin-orbit part. Eq. l|^5|) can 



then be written for this particular case as: 

M li 7- 

E E (y'f'''m\o^\y'^''"^^^ 

7^' J"M" 

J OA|3;^Aio(J^^)) . (31) 

It is important to keep in mind that the operator 
has three terms (see Eq.®), each expressed in one of 
the three sets of Jacobi coordinates. When inserted in 
Eq. ([?Tl) . the two terms in differing from the set of 
Jacobi coordinates i must be transformed into this set 
by use of Eq. © . 

When only one of the three particles is charged the 
operator reduces to one term. A partial wave de- 
pendence in the interaction between the charged particle 
and any of the other two will produce a non-vanishing 
contribution to the energy-weighted sum rule according 
to Eq. pTjl . However, if the only partial wave dependence 
appears in the two-body potential between the two neu- 
tral particles then s[^''^—0. This is because the oper- 
ator then automatically is written in the same Jacobi set 
as the angular functions in Ea. (j3ip . The operator is then 
independent of flx^, and the integral over these angles in 
the last matrix element of ([?T|) vanishes unless ^" . = . , 
■s" — ^x i s-'^d j". = j'^. and therefore the full matrix 
element vanishes. 

The integrals over flx^ and Qy. in the two matrix ele- 
ments that appear in Eg. pip can be calculated analyti- 
cally, because Xi = psinai, and therefore the two-body 
potentials are independent of the angles ilx^ and i^y^- 
The expressions for these two matrix elements are given 
as Eq. (jA6p of the appendix for the particular case of par- 
ticles without spin. 

As an example we consider dipole excitations (A=l) in 
a system of three spin zero particles with equal mass 
TO=4TOAr, where rriN is the nucleon mass, and where 
two particles are neutral and one particle has a charge 
equal to twice the proton charge. We consider only s 
and p waves in the calculation. The two-body interac- 
tions are taken to be gaussians {Vs.p{r) = Ws^pe~^ 1^ ) 
with equal range 6 for s and p waves. We have con- 
structed a 0"*" ground state with a very large contribu- 
tion of s-waves, and a few percent of p-waves. This has 
been done by taking 6=2.98 fm, Ws=Wp^-^.\% MeV 
for the interaction between the two neutral particles, 
and Ws = Wy=—\.\% MeV for the interaction between the 
charged particle and one of the neutral ones. The bind- 
ing energy of the 0+ state is —10.38 MeV. According 
to Eq. ([?T|l . since only terms with ^ l'^. contribute 

(4'^=4^=0, J- j;,= 4,, and /t;,,,_=l), the con- 

tribution sf'^^ to the energy-weighted sum rule is pro- 
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FIG. 1: For dipole excitations in a system of three particles 
with equal mass {m—4 times the nucleon mass) , and only one 
of them with a charge (equal to twice the proton charge), the 
figure shows the variation of (Si — S^"^') as a function of the 
strength difference (AVF) between the gaussian s- and p-wave 
two-body potentials of equal range (see text). The solid and 
dashed lines give the results when only the s-wave strength 
and only the p-wave strength are changed, respectively. The 
dotted line shows the same variation as a function of the 
strength difference between the gaussian three-body forces in 
the excited (1~) and ground (0^) states when the three-body 
potentials have equal range, when only the strength of the 
1~ states is changed, and when the two-body potentials are 
.^-independent. 



portional to AW — Ws — Wp, where Wg and Wp refer 
to the strengths of the interactions between the charged 
and the neutral particle. 

In Fig[Tl the dashed line (+ signs) shows S'^'''' as a 
function of AW when the strength of the s-wave poten- 
tial between the charged and neutral particles {Ws) is 
kept fixed and Wp is changed. Since the p-wave contri- 
bution to the ground state wave function is insignificant, 
a small variation in the strength of the p-wave poten- 
tial only slightly modifies the ground state radial wave 
functions F^°'^°{p). Therefore, the behaviour of is 
almost perfectly linear with AW. However, if we modify 
AW by keeping Wp fixed while changing Ws, the radial 
wave functions are much more sensitive to a change in the 
s-wave two-body potential, since the s-waves dominate. 
Therefore S^^'' is not a completely linear function of AW 
as seen by the solid line (-x signs- in the figure). When 
AW=0, the total value of 5*1 for this particular case is 
51=9.90 MeV fm^, which means that the contribution 
from AW j^O can be of comparable size, see FiglT] 



C. J'^ -dependent three-body potentials 

When performing three-body calculations it is quite 
usual to employ effective three-body forces to fine tune 



the energies of the computed states. Very often differ- 
ent three-body forces are used to place the lowest state 
with given angular momentum and parity J'^ at the cor- 
rect energy. This means that these three-body potentials 
usually depend on J'^. In this subsection we investigate 
the additional contribution Sp^'' to the energy-weighted 
sum rule arising from this kind of three-body potentials. 
With these assumptions the three-body potential opera- 
tor can be written as: 



36 



(p)|n(j.)J-M)(n(j.)J'^M| (32) 



where n(,/ir) refers to all the additional quantum numbers 
needed to specify each of the three-body states with total 
angular momentum and parity J'^. Following Eq. (|14p we 
can write: 



{no Jo"" Mo \0',^0',V3b\no J^" Mo) 



(33) 



where we explicitly labeled the initial state by its parity 

TTq. 

Substituting now Eq. p2p into (j33p and inserting the 
unity operator between O^^ and in the last matrix 
element, we finally get: 



(34) 



j7rMn(,/ir) 

X {noJ^"Mo\O^J{V,['^^ - V^f^)\n^.j.)rM) 

that gives the contribution to the energy-weighted sum 
rule from J'^-dependent three-body potentials. This con- 
tribution vanishes when the three-body interactions are 
J'^-independent. 

For the special case in which the ground state has Jq — 
the expression simplifies to: 



(36) 



(35) 

where tt = 7ro(— 1)^ and which, except for the difference 
between the three-body potentials, is similar to Eq.®. 
Therefore, the analytic expression of Sf'''^ for Jq—O is 
given by Eg. p^ . but with an additional factor equal to 
the difference between the three-body potentials inserted 
in the last matrix element. 

In Fig[T]the dotted fine (with circles) shows sf^'' for 

the same system and the same transition as for the sf''''' 
case. We have taken the ^-independent two-body poten- 
tials used in subsection IIV Bl as starting point, meaning 
that S'P''^ = Si — 5"^^^ . The result is shown as a function 
of the strength difference (AVF) between the gaussian ef- 
fective three-body forces used to compute the 1~ excited 
states and the 0"^ ground state. The range of the three- 
body force is the same (6.0 fm) for 0+ and 1^. The 
variation in AW is obtained by changing the strength 
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in the three-body force for the 1~ excited states. Then 
the ground state wave function remains unchanged. As a 
consequence, according to Eg. ([34|) . and as demonstrated 



by the dotted hue in the figure, Si 
on AW . The contribution from s[ 



(3b) 



(36) 



depends hnearly 
can be of compa- 



rable size to the value, 51=9.90 MeV fm^, for angular 
momentum independent potentials. 



V. A REALISTIC CASE: DIPOLE 
EXCITATIONS IN ''HE 

The main properties of the borromean two-neutron 
halo nuclei are well reproduced describing them as three- 
body systems made by an inert core surrounded by two 
neutrons. The characteristic feature of these nuclei is 
their large spatial extension, which is responsible for the 
large values of the breakup cross sections after electro- 
magnetic excitation. This can be easily envisaged from 
Eas. p3)) and ((20)) . which depend directly on the size of 
the system. For this reason, electromagnetic excitations 
of two-neutron halo nuclei have attracted a lot of atten- 
tion, specially dipole excitations, which is the dominating 
multipolarity for such excitations. 

In this section we investigate dipole excitations in ^He 
(a+n+n), which is one of the most prominent exam- 
ples of borromean two-neutron halo nuclei. We com- 
pute the three-body states by use of the hyperspheric 
adiabatic expansion method [l^. The neutron-neutron 
and a-neutron interactions are the ones used for instance 
in The computed bound ground state (0+) has a 

two-neutron separation energy matching the experimen- 
tal value of -0.97 ± 0.04 MeV. This is achieved with 
a gaussian effective three-body force with range 2.9 fm 
and strength —7.55 MeV. The continuum 1~ states have 
been discretized by use of a box boundary condition at 
Pmax =50 fm. 



A. Dipole strength function 

The transition probability B from the ground state to 
one of the box discretized continuum states is given in 
Eq.ilT]). To obtain a smooth distribution from the dis- 
cretized continuum we use the finite energy interval ap- 
proximation to the strength function. 



dE 



AB 

'AE ' 



(36) 



where AE is the size of the given energy interval, and AB 
is the sum of the transition probabilities into the states 
whose energies fall into this interval. These values are 
then plotted as function of the central energy values of 
the intervals. 

The interval should be of a reasonable size, i.e. large 
enough to provide a smooth function but small enough 
not to wash out the desired structure. In practice we 
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FIG. 2: Computed (circles) dipole strength function for ®He. 
The curve through the computed points is obtained by cubic 
interpolation. Gaussian three-body forces are used to obtain 
the ground (0^) and excited (1~) three-body states. The 
range of the gaussian (2.9 fm) is the same for both angular 
momenta. The different curves show the results obtained with 
different values of the strength {W'^ ') of the three-body 
force for the 1~ states. The experimental data (shaded area) 
are from [isll. 



have used bins with centers at 0.3, 0.9, 1.5, 2.3, 3.2, 4.2, 
5.4, 6.7, and 8.2 MeV, and a standard cubic interpolation 
to smooth the curve. 

In Fig. [2] we compare with the measured distribution 
the computed smoothed dipole strength function for 
^He for different gaussian three-body potentials in the 
excited 1~ states. A strongly attractive three-body po- 
tential produces a pronounced low-lying peak which for 
even stronger attraction would turn into a bound state. 
Thus the three-body potential should at least be less at- 

(i~) 



tractive than that corresponding to W\\^ ' = —40 MeV 
(&=2.9 fm). For moderately attractive three-body po- 
tentials we observe an increase from zero at threshold to 
a peak value at around I MeV followed by a relatively 
fast decrease towards zero at higher energies. This is 
consistent with the calculations in Q and [l^. Com- 
pared to the experiment [l^ . the theory overestimates 
the strength at around 1 MeV and consistent with the 
sum rule underestimate the strength at higher energies. 
Apparently a significant three-body repulsion in the 1~ 
channel would approach the experimental data presented 
with relatively large error bars. 



B. Sum rule results 

The second and third columns of table [ll] give the non- 
energy weighted {Sq) and energy- weighted {Si) dipole 
sum rule strengths. The experimental data, available 
from (Tsj . are given in the first and third rows includ- 
ing states of energies below 5 MeV and 10 MeV, re- 
spectively. The corresponding theoretical values are ob- 
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TABLE II: Non-energy weighted (So) and energy weighted 
(Si) dipole sum rule values for ^He. The fourth column gives 
the contribution to Si from the kinetic energy operator {S[^''). 
The last column {Si /So) is an average dipole resonance en- 
ergy. The upper and lower part of the table show the ex- 
perimental and computed sum rule results for states below 
energies of 5 MeV and 10 MeV, respectively. The last row 
gives the converged results including all excitations. The ex- 
perimental data are from [l3|. The So values are given in 
units of e fm and Si and S{ are in units of e fm MeV. 
The average dipole resonance energy Si/ So is given in MeV. 





So Si Sj^^ Si/So 


E* <5 MeV (exper.) 
E* <5 MeV (theor.) 


0.59±0.12 1.9±0.4 - 3.22±0.94 
0.66 1.94 - 2.94 


E* < 10 MeV (exper.) 
E* < 10 MeV (theor.) 


1.2±0.2 6.4±1.3 - 5.3±1.4 
1.01 4.43 - 4.39 


Converged 


1.25 8.26 4.95 6.61 



tained numerically directly from the first row of Eq.Q 
and Eq. (|15p , and they are given in the second and fourth 
rows of the table. 

We can see that the computed results for Sq agree very 
well with the experimental values when the sum over the 
excited states in Eq.® is restricted to energies below 
5 MeV and 10 MeV, respectively. As even higher ener- 
gies are included the value of Sq converges to the result 
given in the last row of the table, which agrees with the 
expected result obtained from Eq. fT^ . The converged 
value is already reached with an energy limit of about 40 
MeV. 

Essentially the same happens for ^i. The computed 
values agree reasonably well with the experimental ones. 
Also the result obtained for energies below 10 MeV is 
still clearly below the converged value, which requires in- 
tegration up to energies at least of about 60 MeV. The 
converged value for Si clearly disagrees with the result 
provided by Eq. (|20p . where only the kinetic energy con- 
tribution is considered. This value is given in the fourth 
column of table |TT1 It is important to note that the com- 
puted results given in the table have been obtained using 
the same effective three-body force for the 0"*" ground 
state and the 1^ excited states. This means that the 

(T) 

difference between the converged Si value and S^ is 
exclusively due to the effect of the ^-dependence of the 
two-body a-neutron potentials (see Eq. (PT|) ). As seen in 
the table, this effect is far from being negligible. 

The last column in table [Til shows the ratio between 
the energy and non-energy weighted sum rules, which is 
interpreted as an average energy of the soft dipole mode. 
The value of 6.6 MeV obtained after reaching conver- 
gence in Sn and 5*1 is consistent with previous results, 
like [l^, [21J. where a value of about 5 MeV also is ob- 
tained. In [l^ a clearly smaller value is given (3.8 MeV), 

very likely because they used S^ instead of the full Si 
in the computation of the ratio. In any case these rather 
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FIG. 3: Computed dipole sum rule strength Si for ^He as 
a function of the maximum energy allowed above threshold. 
Gaussian three-body forces are used to obtain the ground 
(0^) and excited (1~) three-body states. The range of the 
gaussian is the same (2.9 fm) for all the states. The different 
curves show the results obtained with different values of the 
strength {W^l ') of the three-body force for the 1~ states. 
The thick curve is the result with the same three-body force 
for the ground and excited states. The experimental data are 
from uM. 



large variations illustrate how important it is in practice 
to use the correct sum rules in such estimates. 

So far all the calculations of the 5*1 strength have been 
performed with the same effective three-body force for 
the ground state and the 1~ excited states. This force 
was adjusted to fit the experimental two-neutron separa- 
tion energy in the O^-state. However, for the l~-states 
the interaction might be different and an additional con- 
tribution to the Si strength would appear as seen in 
Eq.dSSD. 

In Fig. [3] we show the dipole sum rule strength Si 
for ^He as a function of the maximum energy allowed 
above threshold. We have considered gaussian three- 
body forces with a range of 2.9 fm. For the ground state 
(0+) a strength of —7.55 MeV has been used. When 
the excited states are obtained with the same three-body 
force. Si behaves as shown by the thick solid line in 
the figure. This calculation corresponds to the numbers 
quoted in table [TTl When the strength of the three-body 

force used for the 1^ states is changed, S'^''"' depends 
linearly on the strength difference. 

The thin solid lines in the figure show 5*1 for different 

values of the strength {W^l '') of the three-body force 
for the 1^-states. This strength has been changed from 
up to —40 MeV, which is at the limit of producing a 
low-lying narrow 1^-resonance in ^He (understood as a 
pole of the 5-matrix). As seen in the figure, the deeper 
the three-body potential, the smaller the value of Si. 
The converged value can change significantly with the 
three-body force. The result obtained with a strength 
of —40 MeV is about 33% smaller than obtained for 
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wl^l ''=0. Therefore, the value of the average dipolc res- 
onance energy (the ratio between the values of the two 
sum rules) also changes substantially with the three-body 

force, ranging between 7.2 MeV when W^l ^ = 0, and 4.7 

MeV when W^l~^ = -40 MeV. 

As seen in the figure, all the curves agree well with the 
experimental value obtained for states below a maximum 
energy of 5 MeV. When this maximum energy is 10 MeV, 
all the computed curves are below the experiment. Al- 
though this experimental value has a rather large error 
bar, it is clear from the figure that the smaller strengths 
in the three-body force for the 1~ states are closer to 
the experimental value. This tendency is consistent with 
the fact that a 1 "-resonance has not been found exper- 
imentally since a strength weaker than about —30 MeV 
also excludes such a state in computations. Extension 
to include higher energies in the experiment would al- 
low distinction between the values obtained for differ- 
ent strengths. Probably ''He is a very favorable system 
for this investigation because the core excitations are ex- 
pected to be negligible. 



VI. SUMMARY AND CONCLUSIONS 

We have derived general expressions for the non-energy 
and energy-weighted cluster sum rules for excitations of 
three-body systems. We consider transitions arising due 
to electric multipole operators of any order and each 
of the constituent particles (clusters) may have a finite 
charge. The most obvious nuclear applications are in 
systems close to three-body thresholds where three-body 
clusterization frequently seems to be a dominating part of 
the structure. This is also the region where the spatially 
extended and weakly bound halos appear. 

Accurate calculations of three-body wave functions re- 
quire in general decomposition into Faddeev components 
either by directly solving the Faddeev equations or by a 
variational procedure including similar components ex- 
pressed in the different Jacobi coordinates. Derivation 
of the non-energy weighted sum rule only relies on the 
use of a complete set of intermediate wave hmctions. 
Therefore only matrix elements of the multipole oper- 
ators enter into the expressions whereas the interactions 
disappear altogether, except of course indirectly through 
the properties of the excited continuum states. However, 
the properties of the wave functions are essential and in 
particular the different Faddeev components give rise to 
crucial interference effects when more than one particle 
carry a charge. 

Such interference effects are also crucial for the energy- 
weighted sum rule where in addition also the properties 
of the interactions arc essential. This sum rule is tra- 
ditionally derived as a double commutator between the 
hamiltonian and two multipole operators. Usually then 
only the second order derivatives from the kinetic en- 
ergy operator contribute while the potentials including 



the spin-orbit terms commute with the multipole opera- 
tors and lead to vanishing contributions. However, when 
the interactions are angular momentum dependent, the 
double commutator does not vanish because the multi- 
pole operators themselves also carry angular momentum. 
These contributions must therefore be computed and in- 
cluded in the sum rule estimates. Still the character of 
sum rule remains in the sense that no matter how the ex- 
citations arc distributed, they must add up to the value 
given by the sum rule which only depends on properties 
of the ground state and the interactions. 

The angular momentum dependence and the subse- 
quent contributions to the sum rule are separated into 
terms arising from the two- and three-body potentials 
which often in accurate three-body computations depend 
on angular momentum. A possible sequence to determine 
appropriate potentials could be first to adjust the two- 
body potentials independently for each partial wave to 
known two-body bound or continuum properties. Second 
to fine-tune the three-body state computed with the two- 
body potentials to a desired energy by adding a short- 
range three-body potential with as little structure as pos- 
sible in order to maintain the properties provided by the 
two-body interactions. Both types of angular momen- 
tum dependence are important as they turn out to give 
substantial contributions to the ordinary kinetic energy 
contribution to the sum rule. 

To assess numerically the relative importance of these 
new sum rule contributions we investigate the electric 
dipole excitations of the ground state of the well known 
halo nucleus ^He. We first notice that the strength dis- 
tribution has a peak at around 1 MeV, falls off at higher 
energies and in practice reaches zero at about 60 MeV. 
The contribution to the energy weighted sum rule from 
the two-body potentials amounts to 2/3 of the kinetic 
energy contribution. 

The contribution from the three-body potential de- 
pends on an expectation value of the difference between 
those potentials for ground and excited state angular mo- 
menta of and 1, respectively. Thus for state indepen- 
dent but finite three-body potentials we arrive at the es- 
tablished result of zero contribution. However, the sen- 
sitivity to the difference in these three-body potentials 
is significant. Realistic potentials give estimates of up 
to 30% of the kinetic energy value. This is then also an 
estimate of the sensitivity of the soft dipole mode to the 
three-body potential. 

In conclusion, we have generalized the energy weighted 
and non-energy weighted cluster sum rules for electric 
multipole transitions to angular momentum two- and 
three-body interactions. The additional contributions 
can be comparable in size to the ordinary terms arising 
from the kinetic energy operator. 
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APPENDIX A: INTERMEDIATE EXPRESSIONS 
AND FORMULAS 

In this appendix we give some of the intermediate ex- 
pressions obtained when deriving Eq. (jl9p . 

The commutator between the electric multipole opera- 
tor O^t and the kinetic energy operator T can be written 
as: 



Writing now the two spherical harmonics in terms of a 
single one, and summing up three of the 3-j symbols [lB| 
one gets: 



[[0',\T],0',] = ^J2 (c,)'(cO'^.^feEE^/2ATT (A5) 

— 1 i—l m,n 



£ \-£ A \ / (2A + 1)! 

m n -1^ ) \j i2£ + l)\{2X-2£ + iy. 

(sin 93,fc)^-'^(cos ifi,kyV£{2£ + 1)VX{2X + l)yt^''^ 
^(_l)i-M+--n^(2A - 1){2£ - 1)(2A + l)Yl,,{ny,) 



\-l £-1 A \ \ £ A\l A £ A 
0/l^-m-ni1^-lA-ll 



[0',\T]^-Y,Zk^[\rk-R\'Yl^{nyJ,Ay,] 



k=l 



= E — ^y>^ (l^*^- - R\^YxA%J) ■ V,, , (Ai) 

where we have used that Ak{rj.Yj,jn{^k))=0- 

The scalar product in Eq. fTSl) can be made by use of 
the gradient formula. A derivation of this formula can be 
found for instance in chapter 5 of [22i |. from which one 
has 



v(0(r)y,„(r!,)) 



21+1 I \dr r 



1/2 



2£+lJ \dr r 



d t+l 



(r)y^,f_i^,„(r!),(A2) 



where 



yj-,,™(i7) = E^^^'"(^)(^'"^;i''^l^''"^)^9 (A3) 



with eo ~ 62, and e±i — =F(e2; ± ie.y)/\/2. 

Use of this expression permits to rewrite Eq. (fT8|) as: 



[[0^t,T],0^] = E (cO'(cfe)'2>^fcEE^^2ATT 



X~£ A 



47r(2A + 1)! 



m n J \l (2£+l)!(2A-2^+l)! 

a;^~'^FA-£,n(!^;cfc)(sini^ifc)^"'^(cosvp,fc)'jy^+*~^ (A4) 
(2A + 1){2£ + l)VAE(-l)'n-i,,(t^,JF.-i,.(n,J 



(-1) 



X- 1 1 X \ £- 1 1 



which after summation over m and /i (the summation 
over /Lt comes from ((M])), leads to Eg. p^ for the con- 
tribution of the kinetic energy operator to the energy- 
weighted sum rule S\ 

We close the appendix giving analytical expressions 
for the two matrix elements entering in Eq. (j3ip for the 
particular case of particles without spin. The expressions 
are obtained performing analytically the integrals over 
and Oj,; : 

(j;f ^0^0 (no |o;^+ \y%'}':'-'" m) = (-i)^+^o+Mo+4', -I'l^ 



(2L" + l)(2Lo + l)(2g^ + 1)(24, + 1){2£'^^ + l){2£y^ + 1) 



47r 



L" Lo X \ ^r* 



-A/" Mo M y " ^Z^V 2^K2A-2^)! 



2A! 



/ V 



L" Lo X 

£y^ ^yi ^ 



COS ai I 



fc=l -^0 
P^„ * (cos2Qi)"i^i t 2 2 (cos2ai) 

{smipki) {cos ipki) (sma,) (cosaO , 



(A6) 



where N ^ " is the normalization constant of the hyper- 
spherical harmonic yf^f^^Q), whose precise form can be 
found for instance in [l3l |. 

The expression for the second matrix element in 
Eq. ([3T|) is identical to Eq. (jA6p but with the function 

""''(psinai) — ^"''(psinai) ) included as a factor 



in the integrand, and with primes on the quantum num- 
bers K, l^. and ly.. 
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